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values of the heliocentric distance. The model and early results 
are discussed herein. 
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Prograa of Besearch in the Third Quarter 


During this past quarter, research activities have concentrated on 
development of the algorithm for Including collisions in the Monte Carlo 
particle-trajectory model. 

Multiple Collisions in Cometary Atmospheres 

Models of neutral cometary comae have generally dealt with the two 
extreme cases: one where the densities are high enough so that all neutral 

species are thermalized and are characterized by a bulk radial flow away from 
the nucleus, and the other where the densities are so low that essentially no 
collisions occur. A characteristic distance from the nucleus separating these 
two regions has been defined where the mean free path for a collision is equal 
to the distance from the nucleus. Whipple and Huebner (1976), for example, 
express this collision radius as 


r -i 22 

c 4irv 

where Q * molecular production rate 

— 15 2 

o - cross-section for elastic scattering ~10 cm 
v « molecular outflow velocity. 

The value of r c varies from ~10^ km for a bright comet like Kohoutek near 
perihelion, down to ~10^ km for a small short period comet like P/Encke at 
1 AU. 

In reality, though, there is no clear separation between the two zones. 
In order for a typical radical (CN, C2, OH) which has been emitted from its 
parent with some exothermic velocity to be truly thermalized, it needs on the 
order of 10 collisions with outflowing water molecules. On the other hand, 
many radicals produced outside of r £ will be subjected to at least 1 to 5 
collisions, thus questioning the true applicability of a free flow model. 

What we truly have, then, is a large regime on the order of 0.2 r c to 5 r c 
where partial thermalization will occur. 

The true random walk nature of the elastic scattering problem is modeled 
in a straightforward manner in the Monte Carlo particle-trajectory model 
(MCPTM). Kitamura et al. (1983) have included collisions in a Monte Carlo 
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model essentially based on our original model (Combi and Delsemme 1980a ). 
Since they followed the correct procedure, we have adopted a similar 
algorithm. There is one problem with their procedure which involves the 
expression they assume for the density of outflowing molecules with which the 
modeled neutral radicals are colliding. They assume a Haser-like description 
for the density which is 


n H 2 0 


(r) 


Q 1 
Airv r 7 


ex P (-r/vT H20 ) 


where Q, v are as above 

r * distance from the nucleus 
t^O ■ photochemical lifetime for H 2 O. 


Strictly speaking, this is the correct density distribution for H 2 O molecules, 
but as the t^O molecules dissociate, they produce OH and H, which still 
provide targets with which radicals collide. Since the effective cross- 
sectional area for the fragments is larger than the parent ^0 but their 
outflow velocity is, on the average, larger, a better approximation would be 
to simply drop the exponential term in the density expression. This yields 


"H 2 0 (r) * W* . 

For the simple case of no radiation pressure, this yields a closed form 
algebraic expression for the collision length of an arbitrarily moving neutral 
species, in place of the integral found by Kitamura et al. 

The random walk algorithm begins with the production of a neutral radical 
by its parent with a given velocity vector. From this location and velocity 
vector, a collision length is given by assumed values for the collision cross 
section, the molecular (^0) production rate, and a random number. The 
neutral is then displaced this distance, scattered elastically off a moving 
water molecule and continued on for another collision. The trajectory 
proceeds this way until the observation time snap-shot. 

Preliminary model runs have been performed for the case of CN produced by 
the photodissociation of HCN (Combi and Delsemme 1980b). We have chosen 
parameters appropriate for the photochemical lifetimes at 1 AU and have 
neglected radiation presuure, which is of less importance at 1 AU, for the 
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purpose of simplifying the preliminary analysis* Subsequent models will test 
both processes simultaneously* The physical parameters of the coma models are 
listed in Table 1* 

Models were run for cases of general molecular production rates from 
1 x 10 2 ® s~ l to 1 x 10 30 s' 1 . The collision statistics from these model runs 
are summarized in Table 2. As can be seen, for larger gas production rates, 
only half of the CN radicals produced are in a true free-flow regime. 
Furthermore, the effect collisions have on the determination of radical 
production rates is also evident in these model runs. For the case of the 
largest gas production rate, if one neglects collisions in the model, an 
underestimate of 11% to 15% in the determination of radical production rate 
would be made, depending upon the size of the photometer aperture used. The 
systematic errors, when reducing data for even smaller heliocentric distances 
(r fl ~ 0.5 AU), should be quite important in bright comets such as P/Halley, 
since at these distances the collislonal effects are compounded by larger gas 
production rates and smaller radical parent scale lengths. 



Program of Research for the Fourth Quarter 


Research efforts in the last quarter will concentrate on (1) model 
analysis of the effect of collisions on the spatial distributions of the 
neutral radicals, ( 2 ) evaluation of newly published brightness profiles of 



C 2 (Cochran 1985), (3) incorporation cf time dependent production rates and 
lifetimes into the model code, and (4) preliminary model runs for OH, C 2 
and C 3 . 
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Table 1 


CN Radical Coma Model Parameters 


r H " 

1 

AU 

t(HCN) - 

9 

x 10 4 s 

x(CN) - 

2. 

75 x 10 5 

v( parent) ■ 

0 . 

58 km/ s 

v e (CN) - 

1 . 

02 km/ 8 

scattering 

1 

x 10" 15 





